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ABSTRACT 

The problem of an ultrasonic reflectivity tomography is similar to that of a spotlight-mode aircraft 
Synthetic Aperture Radar (SAR) system. The analysis for a circular path spotlight mode SAR in this paper 
leads to the insight of the system characteristics. It indicates that such a system when operated in a wide 
bandwidth is capable of achieving the ultimate resolution; one quarter of the wavelength of the carrier 
frequency. An efficient processing algorithm based on the exact two dimensional spectrum is presented. The 
results of simulation indicate that the impulse responses meet the predicted resolution performance. 
Compared to an algorithm previously developed for the ultrasonic reflectivity tomography, the throughput 
rate of this algorithm is about ten times higher. 

INTRODUCTION 

The geometry of a two dimensional reflectivity tomography can be depicted in Figure l.a. It consists 
of a water tank of dimensions greater than that of the object to be imaged. A moving transducer and a 
receiver are scanning over the outer circle. Ultrasonic pulses are emitted from the transducer. Signals reflected 
from the object are sensed by the receiver and collected for further processing. For simplification, we shall 
consider that the transducer and receiver are collocated. Furthermore, the object shall be an idealized 
reflecting medium in which the velocity of the sound is constant, the medium is weakly reflecting, and 
absorption is uniform over the region of interest [1]. 

A spotlight mode SAR in circular path as shown in Figure l.b is analogous to the reflectivity 
tomography system addressed here. In this SAR system, the altitude is assumed to be much smaller than the 
radius of the flight circle. The pulses transmitted periodically from the radar is a wide band electro-magnetic 
wave traveling at a constant speed in the atmosphere. The signal reflected from the ground spot within the 
radar footprint is received by the same radar used for emitting radar pulses. Echo signals are demodulated by 
the carrier frequency, converted into discrete digital format, and stored in an on-board recorder. The stored 
data will be processed by an on-board or ground based processor to reveal the ground image. The design of 
the radar pulse and echo timing for both systems must be coherent. 
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... Since the analysis and processing algorithm description for such systems were originated from the 
T S , AR SyStem ' U wiU 3150 ** 8 iven in following paragraph. However, it should be noted 
la processin S al g° rithm are applicable to the ultrasonic reflectivity system as shown 


ANALYSIS 

A wide azimu * beam SAR can offer higher resolution or wider azimuth viewing angle; two factors 
at help better characterize the backscattering property of targets for various scientific applications. One 
disadvantage of a wide beam SAR is that a much higher pulse repetition frequency (PRF) is usually required 
since PRF is proportional to the radar beam angle. This problem can be resolved using a spotlight-mode 
concept: steering a narrow beam SAR to a fixed spot on the ground. The drawback of a spotlight-mode SAR is 
its limited coverage. 

A spotlight-mode SAR can be operated from a circular flight path as shown in Figure l.b. This type of 

» t m0d f SAR offers several advantages: (1) relatively easy to achieve the ultimate azimuth resolution 
(2) it allows a full 360 degree of viewing angle, (3) there is no need to steer the radar, and (4) the required PRF 
can be scaled down according to the ratio between the radius of the radar spot to the radius of the flight path. 

k For an aircraft SAR in a circular path as shown in Figure 2, the slant range history of a point-target is 

R(9 i) = VA + R-cos0i 

where 

A = Rq + 2 R b - 2RoR h cos 8 e , and B = 2R b {RQ cos 6 e - R b ). 

In the above equations, 6 e is the radar elevation angle, Rq is slant range between the radar and the target at 
the beam center, R b is the radius of the flight path, and 8 X is the angle between the radar at the minimum 
range to the target, the center of the flight circle, and the radar at time of interest. Since time is directly 
proportional to 9 lf this function can also be expressed using time as the variable. 


R(t) = A + B ■ cos(— t) 

V R b 

where v is the speed of the aircraft. Let Z denotes file altitude of the aircraft. It is obvious that Rq and Z are 
related by Z = Rq sin 9 e . Z=0 is a special case that exactly follows the geometry of the ultrasonic reflectivity 
tomography as shown in Figure l.a. Below, we shall denote the distance between the target and the center of 
the path projection on the ground as R T , which is equal to R b -Ro cos 9 e . 

Minimum 



139 



To provide better insight into the spot-light SAR characteristics, it is assumed that the radar beam 
width is unlimited such that any target within the path projection is illuminated all the time during mapping. 
To determine the azimuth resolution, it is necessary to determine the maximum Doppler bandwidth. This can 
be accomplished by evaluating the zero crossing time of the Doppler frequency rate. The Doppler history is 
the derivative of the slant range, i.e. 




-d(2R(t)/X) 

dt 




where *o is the wave number and is equal to 2n/X. X is the wave length corresponding to the center 
frequency of the radar pulse. Plots of both the slant range history and Doppler history are given in Figure 3a 
and b. The history of the Doppler frequency begins with OH z at the minimum slant range point, decreases to 
the minimum Doppler (negative value), increases back to 0Hz at the maximum slant range point, keeps 
increasing up to the maximum Doppler, and then decreases back to 0Hz after a complete cycle. Based on this, 
one can divide a complete circle into two apertures with equivalent bandwidth. The aperture with a Doppler 
ranging from its minimum to its maximum and consisting of the minimum slant range point is referred to as 
the principle aperture. The other aperture spanning the rest of the flight path is referred to as the 
complement aperture. 

V»> 



The Doppler frequency rate variation is given by 
The solution of the time at which f r equals zero is given by 



^cos-'r^- 82)1,2 > 
v B 


Therefore, the Doppler bandwidth, denoted as F, is given by two times the absolute value of the Doppler at 
one of the solution given above. Since F varies as a function of both the aircraft altitude Z and target radius 
R T , it is given by 


F(Z,R t ) = ?j^-(A 2 - - B? - A Ja 2 - B 2 I 

The azimuth resolution is thus given by the reciprocal of the bandwidth multiplied with the effective velocity 
(R T /R b )v,or 


AX(Z,Rj-) = (~f-v) 


1 


R b F(Z,R T ) 
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^variations of the resolution and die Doppler bandwidth as a function of the target location are 
plotted in Figure 4.a and Figure 4.b. It should be noted that these plots also include the targets located 
outside of the vertical cylinder containing the orbit because they can also be imaged as long as range 
ambiguity can be avoided. It is interesting to see that the ultimate resolution of A/4 can be achieved for 
targets falling on the orbit plane and bounded by the circular orbit The resolution width increases linearly as 
targets shift outside of the orbit in the radial dimension. The Doppler bandwidth decreases as the targets 
approach to the center of the orbit. This implies that for targets located cm the orbit plane and bounded by the 
orbit, the required PRF is proportional to the distance from target to orbit center. 




Figure 4.a Resolution vs target location Figure 4.b Bandwidtn vs target location 

Consider a special case where Z = 0 and R T approaches R b . According to the bandwidth equation 

given above, we may find that the required PRF is equal to 4v/A . This indicates that the required sampling 

spacing is exactly A/4 which equals to the ultimate azimuth resolution. The time interval of the principle 
aperture is given by r r 


T (Z,R T ) = ^cos~ l ( 


-1 ,-A + (A 2 -B 2 ) 1 ' 2 


B 


) 


It will be more convenient to express the aperture time interval as a value normalized by the period of a 
complete flight circle. The complexity of SAR processing is usually determined by the number of samples 
within the aperture or the value of the time-bandwidth product TBP given by TBP(Z, R T )=T(Z, R T )F( Z, R T ) 
The time interval and time-bandwidth product of the principle aperture are shown in Figure 5a and 5.b. 


In summary, the analysis in this section leads to the determination of the PRF for operating the 
spotlight radar, the predicted resolution, the time interval of the aperture and the time-bandwidth product to 
be selected for signal processing. r 



Figure 5.a Aperture Interval vs target location 



Figure 5.b Time-bandwidth product vs target location 
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PROCESSING ALGORITHM 


Several algorithms were presented for processing spotlight-mode SAR data or the rotating object 
data. These algorithms include the well known range-Doppler algorithm applied to sub apertures, the 
backprojection processing method (Munson, et. al., [2]) commonly used in computer-aided tomography 
(CAT), and the polar format processing algorithm (Ausherman, et al [3]). One essential assumption required 
for both the CAT and polar algorithms is that the dimension of the imaged area is much less than die radar to 
target distance. Other processing algorithms devised for imaging rotating object include a tomographic 
extension of Doppler processing algorithm (Mensa, et. al. [4]) and a range-Doppler processing algorithm 
(Walker, J. [5]). The first method is suitable only for imaging sparse arrays of objects due to its higher 
integrated sidelobe ratio (ISLR); and both algorithms also rely on the assumption of large radar to target 
distance. An exact solution for a circular aperture acoustic imaging system was presented by Norton [1]. This 
algorithm requires the implementation of a quasi fast Hankie transform which is not very efficient. 

The difficulty in processing spotlight-mode SAR data collected from a circular path is that the 
Doppler history of a point target involves many higher order terms and that the depth of focus is very 
shallow. To overcome the first problem, the algorithm proposed here make use of an exact 2-D spectrum of a 
point target (Jin, 1992 [6]) in a range Doppler like processing approach. To overcome the second problem, this 
algorithm updates its reference function as frequently as required. 


According to [6], the magnitude and the phase of a reference spectrum are given by 


A(co r , ( 0 ) = V2»r 


d 2 (-(<o r + <p 0 )-2R(t,)/c» 
d f 2 


- 1/2 


yr(fi) r ,fi)) = exp{ y( l — ^ ^ - 2 k r jA + Ba(o) r , a) - -j) } 

v/ Rb 4 

where 0) is the azimuth angular frequency, (o r is the range angular frequency, (Oq is the angular frequency 
of the carrier, K r is equal to <» r /c, where c is the speed of light, and 


CD 2 Iffl 4 
a(ffl r .o)=(-^-±J-^ 


/ R b 




it 4 


The energy distribution of the 2-D spectrum of a point target is given in Figure 6. It can be seen that 
the amount of area with energy is proportional to the distance from the target to die center point. 


Range frequency 
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The basic correlation steps for the spotlight-mode processing are summarized below. 

(1) ^l° r ppp b ?r Pling ( ° rthe **? pukeS accordin S to ^ of the spot area. This step is not required 
if the PRF of the sensor is tuned to that radius such that there is no redundant data. Transform the SAR 
data into its spectrum by a 2-D FFT process. 

(2) Perform a 2-D SAR correlation. For each azimuth line with a constant range position, a reference spectrum 
is generated according to equation (1). Correlation is performed by multiplying the data spectrum and the 

conjugated reference, averaging in range, and performing an inverse azimuth FFT. This is repeated for 
each azimuth line. r 


(3) Perform geometric resampling to correct for the geometry and grid spacing. The mapping between image 
pixels generated from step (2) and points on the ground for a spotlight-mode SAR is given in Figure 7. 


Azimuth 


t 

Range 


Image 
Pixel Grid 



Figure 7. Image pixel to projection mapping in a spotlight mode SAR 


COMPUTATION EFFICIENCY 

„ , e'gorithm described above takes the advantage of computation efficiency of FFT operations, 

therefore, it is a highly efficient algorithm as compared to a previous algorithm proposed by Norton [11. To 
illustrate this, the Norton's algorithm is briefly reviewed first. 

Let f(r,0) denotes the reflectivity of the medium and g(r,0) denotes the echo signals received in a 
polar coordinate with radial position r and angle 6. f(r, 9) can be recovered by the following three steps: 

(1) Perform spherical harmonic expansion for g(r, 9) for its coefficients, £ n (r) = — fg( r , 9) e ~ in6 d9 

2n Jo ' 

(2) Get die spherical harmonic expansion coefficients /„(r) of f(r,9) by f n (r)=H„{ 


j is i 8n (P) i i 
J„(R z) ° 2 np }z)r 


(3) Reconstruct f(r,9) from /„(r) by f(r,9)= £/„(r) e 


,in9 


< . u l T, the S6C ? nd Ste Pf Hake's transform of the 0-th order and n-th order are implemented by a quasi 
ast Hankie transform. It is assumed that in both algorithms the reference functions are generated before 
processing and permanently stored in the computer memory. The computation efficiency will be examined 
on y on the amount of complex multiplys performed to yield an image of a diameter of 2 N pixels. These 
results are summarized in Table I and Table II. 
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Table I. Amount of complex multiply in spotlight mode algorithm 


PROCESSING STEP 

2-D FFT FOR 
ECHO DATA 

REFERENCE 

MULTIPLY 

1-D INVERSE 
FFT 

NUMBER OF MULTIPLY 

7tN 2 Log 2 nN 
+7tN 2 Log 2 N 

ItN 3 12 

7CN 2 Log 2 JtN 


Table II. Amount of complex multiply in Norton's algorithm 


PROCESSING STEP 

HARMONIC 

EXPANSIONS 

REFERENCE 

MULTIPLY 

1-D HANKLE 
TRANSFORMS 

NUMBER OF MULTIPLY 

2nN 2 Log 2 nN 

nNN h 

Nk=K\N Ln(K 2 N) 

2 N h Log 2 N h , 1 

N h =K x N L,(K 2 N ) 


The disadvantage of the quasi fast Hankie transform is that data samples must be oversampled from 
N to N h = K X N Ln(K 2 N), where Kj and K 2 are constants. This oversampling increase the amount of data and 
die amount of multiplies. The computation bottleneck of the spotlight algorithm is in the reference multiply 
which is proportional to N 3 . For N being 128 to 512 points, the amount of multiplies in die spotlight algorithm 
is from 1/10 to 1/4 times of that in Norton’s algorithm. Since a large portion of the point target spectrum has 
no energy as shown in Figure 6, the amount of multiply can further be reduced by a factor of two in reference 
multiply and inverse FFT processes. 


SIMULATION 


A simulation was performed to test the proposed algorithm and to verify the analysis on the azimuth 
resolution. Plotted in Figure 7a and 7.b. are the impulse responses of two aircraft spotlight systems, both with 
zero altitude. The first one has a range bandwidth being twice the carrier frequency such that the range 
resolution is comparable to the azimuth resolution. The second system has a range bandwidth of .15 times the 
carrier frequency. The 3-dB resolutions in range and azimuth of the first system are very close to A/4. The 3- 
dB resolutions of die second system are also close to A/4 . However, it has a much worse integrated sidelobe 
ratio. 


Range Bandwidth = 2 x Radar Carrier Frequency 

Full Resolution Processing over the Principle Aperture 




Range Bandwidth = 0.1 5 x Radar Carrier Frequency 


Full Resolution Processing over the Principle Aperture 



Figure 7. a Spotlight SAR impluse response Figure 7.b Spotlight SAR impluse response 
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CONCLUSION 


The spotlight-mode aircraft SAR in a circular flight path has been analyzed in terms of resolution, 
bandwidth, aperture arc, and the time-bandwidth product. An exact and efficient processing algorithm is 
proposed to recover spotlight-mode SAR images. The simulation result of the point-target response indicates 
that for full azimuth resolution processing both the azimuth and range 3-dB resolution can achieve the 
ultimate resolution of X/4, regardless of the range bandwidth. This algorithm is directly applicable to 
ultrasonic reflectivity tomography system. With further optimization in computation efficiency, the 
processing throughput rate shall be about ten times faster than a previously proposed algorithm. 
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